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Abstract 

We investigate nonmagnetic impurity effects in two-band superconductors, 
focusing on the effects of interband scatterings. Within the Born approxima- 
tion, it is known that interband scatterings mix order parameters in the two 
bands. In particular, only one averaged energy gap appears in the excitation 
spectrum in the dirty limit. [G. Gusman: J. Phys. Chem. Solids 28 (1967) 
2327.] In this paper, we take into account the interband scattering within the 
t-matrix approximation beyond the Born approximation in the previous work. 
We show that, although the interband scattering is responsible for the mixing 
effect, this effect becomes weak when the interband scattering becomes very 
strong. In the strong interband scattering limit, a two-gap structure corre- 
sponding to two order parameters recovers in the superconducting density 
of states. We also show that a bound state appears around a nonmagnetic 
impurity depending on the phase of interband scattering potential. 
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I. INTRODUCTION 



Since the discovery of superconductivity in MgB2, [1] great efforts have been done on 
this material in order to clarify the character of this 40[K]-superconductor. Experiments 
on the specific heat [2] and penetration depth [3] indicate a finite energy gap below the 
superconducting transition temperature T^. The presence of a Boron- isotope effect [4] implies 
that phonon plays a crucial role in the pairing interaction. However, the origin of this 
superconductivity is still in controversial and various mechanism have been proposed. [5-7] 

Besides the mechanism, possibility of multi-gap superconductivity has been discussed in 
this material. [8-11] Indeed, a two-gap structure has been directly observed by tunnelling 
experiments. [12,13] The temperature dependence of the specific heat also indicates the 
presence of two energy gaps. [2] Since MgB2 has two distinct Fermi surfaces originating 
from (J- and 7r-orbitals of boron atoms, [14] it is probable that the superconducting order 
parameters are different in the two bands. We also mention that multi-gap superconductivity 
has been investigated also in the superconducting state of transition metals. [15-23] 

When the order parameters are different among bands, it is known that impurities with 
interband scatterings strongly affect the superconducting state even if impurities are non- 
magnetic. [19-23] In particular, Gusman [19] showed within the Born approximation that 
the order parameters in a two-band superconductor are mixed with each other by interband 
scatterings. In the dirty limit, only one averaged energy gap appears in the superconducting 
density of states. Although this impurity effect was originally considered for the supercon- 
ductivity in transition metals, this effect can be expected also in MgB2 if this superconductor 
really has two distinct order parameters in the two bands. 

So far, the mixing of order parameters by interband scatterings has been discussed within 
the Born approximation. [19-23] With this regard, we mention that higher order scattering 
processes beyond the Born approximation sometimes lead to qualitatively new effects in 
the case when the impurity scattering affects superconductivity. (For example, formation 
of a bound state around a magnetic impurity in s-wave superconductivity. [24]) Thus it 
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is expected also in the present case that higher order scattering processes may cause new 
phenomena which cannot be described within the Born approximation. 

In this paper we investigate nonmagnetic impurity effects in a model two-band supercon- 
ductor with two different order parameters. We focus on the mixing of the order parameters 
by interband scatterings. In contrast to the previous works which employ the Born approxi- 
mation, we take into account multi-scattering processes in terms of the inter- and intra-band 
scatterings using the i-matrix approximation. We clarify how the mixing effect of the order 
parameters which has been discussed within the Born approximation is modified by higher 
order scattering processes. We show that, although the interband scatterings are crucial for 
the mixing effect, at the same time, this effect is suppressed when the interband scattering 
is very strong; this result is not obtained within the Born approximation. As another new 
result by the t-matrix approximation, we also show that a bound state is formed around a 
nonmagnetic impurity depending on the detailed character of impurity scattering. 

This paper is organized as follows: After explaining our model in §2, we consider a single 
impurity problem and discuss a bound state in §3. In §4, we investigate the case of finite 
impurity concentration, which is followed by summary in §5. 

II. MODEL TWO-BAND SUPERCONDUCTOR 
A. Hamiltonian 

We consider a three-dimensional two-band s-wave superconductor. We refer to the two 
bands as a- and &-band. The Hamiltonian is given by [15-19] 

H = Y1 ^A^^P'y + ^ViX<j - 9aa Yl ^T^^-Pi^-p'i Vt - 9bb Yl ^it^-Pi^-p'i Vt 

Pi"' Pif p,p' p,p' 

- 9ab Y (^ll(^-pl^-p'lK'l - 9ba Y ^It^-pi^-P'i^P'T 
P,p' P,p' 

+ i^imp. (2.1) 

Here, (bpa) represents the annihilation operator of an electron in the a-band (6-band). 
The kinetic energies of the a- and 6-band measured from the Fermi level are, respectively, 
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given by and e^. In this paper, we do not take into account the detailed band structure of 
MgB2 and simply assume two isotropic three-dimensional bands. The a-band (6-band) has 
an intraband pairing interaction described by gaa {gbb}- In addition, eq. (2.1) also includes 
an interband pairing interaction described by gab and g^a {— g^},)- We do not specify the 
origin of the pairing interactions and, for simplicity, assume the same cut-off energy cuj^ for 
all the pairing interactions. We briefly mention that different cut-off energies between the 
two bands are crucial in considering the isotope effect as discussed by Kondo. [16] The last 
term in eq. (2.1) represents the nonmagnetic impurity scattering: 

-f^imp = J2 e'^"""''^'^* [VaaO'laO'p'a + V<t + VahdlaK'o + V^hl^d^'a] , (2-2) 

where Vom and Vbb, respectively, describe intraband impurity scattering in the a- and 6-band 
while Vab and Vba (= "^afc) represent interband scattering between the two bands. In eq. 
(2.2), we have neglected momentum dependence of scattering potentials for simplicity and 
retained the s-wave scattering component only. In the case of a single impurity problem, we 
put the position of the impurity i = and drop the summation in terms of j in eq. (2.2). 

We employ the mean field approximation and introduce the two superconducting order 
parameters in the a- and 6-band: 

A„ = -gaa Ep(apTa-pi), 

[2.6) 

Afc = -9bbY.p{bp^b_pi). 
Then the mean field Hamiltonian of eq. (2.1) is obtained as (we drop constant terms) 

Hmf = E *i(^>3 - A„)*, + ^ ^ie^Ps - ^b)% + i^imp, (2.4) 
p p 

where = (a^^, a_pj^), and = (ft^p &_pj represent two-component Nambu field operators 

in the a- and 6-band, respectively; pi {i — 1,2,3) are corresponding Pauli matrices. In eq. 

(2.4), Aj {j = a, b) are given by 



Re[Aj]pi - Im[Aj]/,2, (2.6) 



where Aj are related to the order parameters as 

A„ = A„ + ^ A,, 
9bb 

A, = A, + ^Aa. 
The impurity term i^imp under the Nambu representation is given by 

Ri,P,p' 

Here the matrix impurity potentials are defined by Vaa — VaaPzi Vbb — VbbPs and 



^Vab ^ 



Vab 



V 



-V. 



ab 



Re[Vab]p3 + ilm[Vab]l, Vba = ^^Jfo- 



(2.6) 



(2.7) 



(2i 



We can take A^ real and Qab — Qba > with use of a gauge transformation. Complex 
quantities in our model are then A^ = lA^je^'''''' and Vab — I'l'able^^"''- 

B. Clean system 



In this subsection, we briefiy review the theory of two-band superconductivity in the clean 
system [15-18] in order to see how the impurity scatterings modify the superconducting state 
in later sections. The one-particle matrix thermal Green function in the j-band {j — a, b) is 
given by 

1 



^"(,10;^) 



(2.9) 



Here cUm is the Fermion Matsubara frequency. When we calculate the superconducting 

density of states (SDOS) in the j-band, we obtain 



Nj{uj) = ImG'5(, ia;^ ^ a; + i6) = Nj{0)Re[- 



TT 



(2.10) 



where Nj{0) is the density of states at the Fermi level in the normal state. Equation (2.10) 
shows that the energy gap in SDOS is governed by, not the order parameter Aj itself, but 
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Aj defined in eq. (2.6). Since includes both and A^, the superconducting states in 
the two bands couple with each other by the interband pairing interaction gab- 

Substituting eq. (2.9) into the gap equation Aj = gjjTY^p,_^^G'j{,iu!m)i2, we obtain 

r^D A,- /3 



A,- = gjjN^iO) / de ^ tanh ^^s^ + | A,f (j = a,b), (2.11) 



where /? — 1/T is the inverse of temperature. Substituting eq. (2.6) into eq. (2.11), we find 
Ajj — CAa, where C is given by 



1 - 9aaNa{0) /o"° ds ' tanh + |A„|2 

C= V±+\^ ^ . (2.12) 

gaa'^NM de-^^^ tanh l^e- + \K\^ 

Thus when A^ becomes finite, A^ also becomes finite as far as gab 7^ and gbb¥' ^- [15] We 
also find that Af, can be also taken real as well as A^ in the clean system. 
The superconducting transition temperature T^. is determined by [15,17] 

^ 2-fUJB r 9aa + 9bb- J (gaa " gbb^ + 4^9ab9bal , ^ ,^ ^ 

= exp (7 = 1.78), (2.13) 

'- ^[9aa9bb " 9ab9ba) 

where gaa = 9aaNa{0), gbb = QbbNaiO), gab = QabNbiO) and gba = 9baNa{0). 

Temperature dependence of the order parameters in the clean system is shown in Fig. 
1. When (jab ~ {gaa, Qbb) (Fig- !(&)), the overall temperature dependence is similar to the 
temperature dependence of the order parameter in an ordinary single-band BCS supercon- 
ductor. In weak couphng BCS theory, the magnitude of the energy gap at T = is related 
to Tc as 2ABcs(r = 0)/Tc = 3.54. On the other hand, as proved by Soda and Wada, [17] 
2A„(r = 0)/% is larger than 3.54 while 2Ab(r = 0)/% is smaller than 3.54 in Fig. 1(a). 
When gab <S {gaa, gbb) (Fig- 1(b)), the smaller order parameter (Ab) shows a curious tem- 
perature dependence: [15] When gab = 0, the 6-band remains in the normal state down to 
T/Tc = 0.513 (where Tc is the transition temperature in the case of gab = 0.002). Thus A^ 
between 0.513 < T/Tc < 1 in Fig. 1(b) is due to a proximity effect caused by the inter- 
band pairing interaction gab- This induced order parameter appears even if the intraband 
interaction in the 6-band is repulsive {gbb < 0). We shows an example in Fig. 1(c). 
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III. BOUND STATE AROUND A SINGLE IMPURITY IN TWO-BAND 

SUPERCONDUCTOR 



A. Equation of bound state 



(3.2) 



Let us consider the electronic state around an impurity at j = 0. In this case, within 
the neglect of the spatial variation of the order parameter, the one-particle matrix thermal 
Green function in the a-band Ga{,',iuJm) has the form 

Ga{: = G% lUJm)S^,^' + G% ZW^)4„(ia;)G'°(', lUJm). (3.1) 

Here the t-matrix iaa obeys 

iaa{iuJm) = ^aa + ^aa Ep G% ia;^)4a(i'^m) + Vab Ep Gg(, ia;^)4a(i'^m) , 

ibaii^m) = ha + Vba Ep G% iuJm)i aai^^m) ~^ ^bb Ep ^bi'' ^^m)iba{^^m) ■ 

The solution of eq.(3.2) is given by 

4a(i'^m) = [1 - Vaa{iuJm)J2'^aai^'^rn)]~^Vaa{iuJm), (3.3) 

P 

where the effective impurity potential in the a-band defined by Vaai^^-^m) includes interband 
scatterings between the two bands as 

Ka(i'^m) = Vaa + VabJ2'^bi^'^rn)['i- - VbbJ2'^bi^^'^rn)]~^Vba- (3.4) 

p p 

A bound state caused by the impurity scattering is obtained as a pole of the analytic- 
continued t-matrix of eq. (3.3). The equation of the bound state energy uj is then given by 
(a;<Min[|A„|,|A6|]) 

a;2 - cos(2^„b) A^Afc 1 



[1 + + a,, + {\Vab\ -QiaOift)), (3.5) 



where = 7i'-^a(0)t'aa and ai, = T:Nb{0)vbb represent intraband impurity potentials normal- 
ized by the density of states; Vab = tt^^ Na{0)Nb{0)vab — \vab\^^^"-'' ■ Since we consider a single 
impurity problem and neglect spatial variation of the order parameters, we have taken both 
Aa and A^ real in eq. (3.5). The derivation of eq. (3.5) is summarized in Appendix A. 



B. Conditions for the existence of bound state 

Since the present theory is symmetric in terms of the a- and 6-band, we may consider 
only the case of | Aa| > | Abl without loss of generahty. The left hand side (LHS) in eq. (3.5) 
equals — cos(2^a(,)sign[AaAf,] at a; = while it becomes ±00 at a; = lA^I depending on the 
sign of A^ — cos(26'a6)AaA6. We then find the following three cases about the formation of 
the bound state: 

(1) cos(2^ab)sign[AaAf,] < 0: A bound state is obtained if the right hand side (RHS) in 
eq. (3.5) is larger than | cos(26'aft)|. 

(2) cos(2^^ab)sign[AaAft] > and A^ — cos(26'a6)AaA(, > 0: A bound state exists because 
LHS in eq. (3.5) monotonously increases from — cos(2^a;,) to +00 as uj changes from 
to lAfel. 

(3) cos(2^ab)sign[AaA6] > and A^ — cos(2^ab)AaA6 < 0: A bound state is absent below 
|Afe| because LHS in eq. (3.5) is always negative. 

A bound state is not obtained when the interband potential scattering Vab is absent because 
RHS in eq. (3.5) goes infinity. In addition, when A^ and A^ have the same sign and the 
interband potential scattering is real, one does not obtain a bound state. Thus, the phase 
of the interband potential scattering and the difference of the signs of A^ and A;, are crucial 
in obtaining the bound state. 

Although the interband scattering is crucial for the bound state, eq. (3.5) also shows 
that the bound state disappears in the limit \vah\ — 00. (Note that RHS in eq. (3.5) goes 
infinity in this limit.) When we write the interband scattering term as [vah<Aa^O(T + h.c.\ {a\^ 
{boa)'- creation (annihilation) operator of an electron in the a-band (6-band) at = 0) and 
diagonalize this term, we obtain two energy levels E± — ±\vab\- In the limit \vab\ 00, the 
lower level E_ is occupied by two electrons while the higher state E_^. becomes irrelevant 
owing to its extremely high energy. Then since electrons no longer come to the impurity 
site, the effect of the interband scattering disappears in the limit \vab\ 00. 
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Dependence of the bound state energy on \vab\ is shown in Fig. 2. A bound state 
appears when the interband scattering becomes finite; however, the bound state energy 
again approaches the energy gap | A^l when \vab\ ^ 1 because of the reason explained above. 
In Fig. 1(a), the bound state energy increases as the phase of interband impurity scattering 
9ab decreases. In the case of Fig. 1(a) (A^ — 2 and Af, — 1), the bound state disappears 
for 9ab > 7r/6 because the case-(3) discussed above is satisfied. The bound state energy also 
increases when the intraband scattering (ab) becomes strong, as shown in Fig. 1(b). This is 
because electrons cannot come to the impurity site when the intraband impurity potential 
is strong; thus the effect of the interband scattering is weakened by the intraband impurity 
scattering. We also find this effect in eq. (3.5); when — > oo and/or — > oo, RHS in eq. 
(3.5) goes infinity as in the case of \vab\ — 0. 

IV. ORDER PARAMETER MIXING IN DIRTY TWO-BAND 

SUPERCONDUCTORS 

A. Formulation: i-matrix approximation 

In this section, we consider the case of finite impurity concentration using the i-matrix 
approximation. After taking spatial average in terms of the distribution of impurities, we 
obtain the renormalized one-particle matrix thermal Green function in the a-band as 

1 



Ga{, i'^m) — 



i(^m - £pP3 + A„ - E„a(i<^m) 
^ ^. (4.1) 

ia;« - e;p3 + A„ 

The self-energy Tiaaiii^m) within the t-matrix approximation is given by 

^aa{iuJm) = nimp[l - 14a(i'^m) Ep Ga{, i<^m)]"^ Ka, 

(4.2) 

Vaa{iuJm) = Vaa + Vab Ep Gb{, iuJm)['i- - VbbGb{, iuJm)]~^Vba- 

Here nimp is the number density of impurities. In eq. (4.2), Vaa has the same structure as 
Vaa in eq- (3-4) except that the renormalized Green function is now used in place of G^. 
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The renormalized Green function in the 6-band is also given by eqs. (4.1) and (4.2) where 
the indices 'a' and '6' are interchanged. We mention that has the same structure as 
in eq. (3.3). 

In Eq. (4.1), (D^ is a renormalized Matsubara frequency and 



Re[A,]pi - lm[A,]p2 



(4.3) 



with Aj being a renormalized A^. [25] Their expressions are derived in Appendix B and we 
only show the results here: 



Ub r \Vba\ 1< 



Db'l + al Qa 



+ r)b- 



Qb 



(4.4) 



Uar At 



b |'6 

/\b — i\b + ttI -, , .9 ?r + 



where Oj = \J ujir? + | AjP, Uj — niinp/(7rA^j(0)) and 

|4 



Vab 



(1 + a|)2 



+ (ttft 



•^6 I- |2\2 
^j— 2 K^«6 ) > 



1 + aV 



(4.5) 



(4.6) 



= 1 + r/a + 

Db^l + rib + 



1 



1 + fla^b 

2 1 

1 + ai 



(4.7) 



B. Gap equation 



The order parameter in the j-band, A^, is determined by the gap equation A^ 

QjjT Y.p,ujm ( , i'^m) 12 : 
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^^=9nli: I . . U = a,b). (4.8) 

As usual, we introduce a cutoff in calculating the summation of the Matsubara frequency 
in eq. (4.8). 

When the interband scattering Vab is absent, a;^ and satisfy 

(4.9) 

In this case, the effect of impurity scatterings in eq. (4.8) vanishes as in the case of the 
nonmagnetic impurity effect in single-band s-wave superconductivity. [26] As a result, super- 
conductivity is unaffected by impurities. The same situation also occurs when (1) intraband 
scatterings or (2) interband scatterings are very strong: In the former case, since the first 
terms in [• • •] in eqs. (4.4) and (4.5) are absent in the limit aj — >• oo, eq. (4.9) is satisfied. 
Physically, when the intraband impurity potential is strong, electrons are excluded from the 
impurity site. Then the effect of the interband scattering is effectively weakened, so that the 
situation becomes close to the case of Vab = 0. In the case (2), since rjj ~ \vab\* for l^a^l ^ 1, 
the first terms in [• • •] in eqs. (4.4) and (4.5) again vanish in the limit Vab — > oo. As discussed 
in Sec. 3.2, the impurity level becomes ±oo in this limit, so that the interband scattering 
becomes irrelevant. As a result, the situation in the limit Vab — > oo is equal to the case of 
Vab = 0. 

Within the Born approximation, the effect of the interband scattering is controlled by 
the damping rate Tj = nrii^npNj {0)\vab\'^ (j = a, h) when the intraband scattering is absent. 
[19] Then as shown by Gusman [19] only one averaged energy gap appears in SDOS in 
the hmit Vj — > oo (j = a,h). [19] (We reproduce this result based in the present model 
in Appendix C.) Namely, within the Born approximation, the interband scattering simply 
promotes the mixing of the order parameters. On the other hand, the present results show 
that this does not always occur; in the limit Vab ^ og, although — > oo, a two-gap structure 
corresponding to | Aa| and | A^l should be observed in SDOS as in the case of Vab = 0. 

In solving eq. (4.8), we employ the following prescription: When we separate eq. (4.8) 
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into the real and imaginary parts, we obtain four equations. However, since we can always 
take Aa real, the number of variables which we have to determine is three, i.e., A^, Re[Ab] 
and Im[A6]. Thus we have to eliminate the redundancy in eq. (4.8). For this purpose, we 
rewrite the gap equation as follows: Prom (4.4) and (4.5), we find that a;^ and A^ can be 
written as 





- I 



(4.10) 



where A ■ (j = a, 6) are related to and A^ as 



A 



CaaAa + ^CafeAft, 
9aa 



a 
1^ I 



In eq. (4.10), is given by 



In eqs. (4.10) and (4.11), Q and (jj> = a,b) are given by 

Db F>al + ab 

- 1 _ _H6_Jv6aP_ J_ 

Da Dbl + a^ ila 



(4.11) 



(4.12) 



(4.13) 



Caa 
Cab 
Cbb 



Ub 9ba Ua 

'-'ft*- 



1 - [;^^a - 



vlb 



9aa 


1 + al 


(jbb Ua 


<b 1 


9ab Da 




9ab Ub 




gbb Db 


l + a^J 


(jaa Ub 





r - 1 _ 

Substituting eq. (4.10) into the gap equation (4.8), we obtain 

/ -1 _ TI" v-^ / N Caa _ 9ab / \ Cat 

^- 9aa^2^s{uJm)^ -^ag — "5 s(^m)^ 



_ 9ba ^ / \ ^"a -I _ " 

-S-bb— 75Z^"5(cc;^)^ l-9bb-3 2^s 



ba 



71 



9aaP 



bb 



b I 



(4.14) 



(4.15) 
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Here Cl'j = y Cj^'m + l^jl^ s{u!m) = F^/\Fuj\. Thus, the determinant of the 2 x 2-matrix 
in eq. (4.15) must be zero in order to obtain finite Aa and A;,: 

(1 -^aa^I]s((^J^)(l -^M^^s(a;J^) -^„,,^ba^I]s((^j&^ " ^I]^(^m)fe = 0- 

(4.16) 

When eq. (4.16) is satisfied, the "eigen- value" (A^, A^) in eq. (4.15) obeys the equation 

A, ^ _ (4.17) 

A„ _ S-abTT^ Cafe 



9hhP "6 

The real and imaginary parts of Eq. (4.17) give two equations, while eq. (4.16) gives only 
one equation because it is a real equation as proved in Appendix D. Thus we can determine 
the order parameters Aa, Re[Ab] and Im[Ab] unambiguously from eqs. (4.16) and (4.17). 

In the clean system, we can take C,j — — F,^ — \. Then the gap equation (4.15) is 
reduced to eq. (2.11) when, as usual, (1) we transform the summations in eq. (4.15) to 



complex integrations ^, (2) change the variables of the integrations as z — > £ = y 2;^ — | Aj 
and (3) introduce a cutoff loy^. 

C. Transition temperature 

The transition temperature T^. is determined by eq. (4.16) with Aa = A^ = 0: 



- 9a,malY.^>rn)T^^ " ^E^^K)^- (4-18) 

III. 

Here C° and C°j, are given by eqs. (4.13) and (4.14) with Aa = A;, = 0: 



2 



1 r 



or i^°i + «j^K|+r,' 



.0 . r'^a Ub \Vba 



(4.19) 



C = 1 - [T^^a 
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/■O 


= 1 — 




Sob 





[^^Vb- 


r° 

S.66 


= 1 - 


r Ua 


S6a 


= 1 - 





9ba Ua 



Qaa D^^l + al \uj,n\ + Tfc ' 



^^2 



9bb_Ua__Vab_ 

gab + 

S-ab Ub via 



b 



1 



^gg Ub via i 1 

In eq. (4.18), ^^(u;^) = where is given by 



(1-^ 



(4.20) 



Ub 


Vb ^ 


Vab 


UaUb 


1 








+ Tb' {l + al){l + al)DODU\oOm 




+ Tb) 



(4.21) 



In obtaining eqs. (4.19)-(4.21), we have used the result that eq. (4.4) at T^. is reduced to 

|2 



• - a • I • ^« r ^"-^ I I 

i(^„ = iu;^ + i— [— -2 +r/aj 



• -b • , • r ^f"* 







+ ira 












+ irb 









(4.22) 



where 



I — 1 2 

^-1 , Pab \2 , / Q;b I- |2\2 



D 



\Vba\ n2 



l + «b 



/I I l^bal \2 I / 



(^a ,_ |2\2 



(4.23) 



In the present case, effects of impurity scatterings are controlled by (i) the impurity 
potential {vaai Vbb, Vab) and (ii) the impurity concentration described by Uj = nimp/(7riVj (0)). 
In Figs. 3 and 4, we show their effects on T^. 

(i) Figure 3 shows the effects of the intcrband scattering Vab on the transition temperature 
Tf.. When Vab becomes finite, Tf. decreases with increasing \vab\ in the region < 
\vab\ ^ 1; however, again increases when \vab\ ^ 1 reflecting that the interband 
scattering becomes irrelevant in the limit Vab — ^ oo. When the interband pairing 
interaction gab is strong, the superconductong states in the two bands strongly couple 
with each other even in the clean system. In this case, even if transfers of Cooper- 
pairs occur between the two bands by interband scatterings, this effect would be weak 
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when the phase of the interband scattering 6ab is absent. Indeed as shown in Fig. 
3(a), the suppression of by Vab is weaker for larger (jab- However, when the phase 
of Vab is non-zero, the Cooper-pair acquires the phase 2 x 6*06 in the transfer between 
the two bands. Then since the phases of the order parameters in the two bands are 
the same {gbb > 0), the appearance of the Cooper-pair having this extra- phase would 
destroy superconductivity. Thus Tc is more suppressed when 9ab is finite in Fig. 3(b). 
On the other hand, since the intraband scattering weakens the effect of the interband 
scattering, the suppression of becomes weak when the intraband scattering becomes 
finite as shown in Fig. 3(c). 

(ii) Effects of the impurity concentration on is shown in Fig. 4. We find that the 
suppression of Tc continues to exist as the impurity concentration increases, which is 
in contrast to the case of the increase of \vab\ in Fig. 3. In the case of 9ab — (Fig. 
4(a)), the suppression of Tc is weaker when the interband pairing cjab is stronger. On 
the other hand, the opposite tendency is obtained in the case of = 7r/2 as shown 
in Fig. 4(b). In this former case, when Qab is strong, the superconducting states in the 
two bands strongly coupled with each other as noted in (i). Thus the mixing effect of 
two order parameters by interband scatterings is weaker for larger gab in the case of 
9ab — 0. However, in the case of Fig. 4(b) {9ab 7^ 0), since the Cooper-pair acquires the 
phase 29ab in the tunnelling by the interband scattering, the depairing effect becomes 
more remarkable when gab is stronger. (In the extreme case gab — 0, one can freely 
choose the phases of and A^. Then the depairing effect by the tunnelling of the 
Cooper-pair by the interband scattering can be eliminated by choosing the difference 
of the phases of A^ and A^ as 29ab- Indeed, one can eliminate 9ab from the present 
model by an appropriate gauge transformation. This means that the depairing effect 
by 9ab is absent when gab = 0.) As a result, in contrast to the case of 9ab = 0, the 
suppression of Tc is more remarkable in Fig. 4(b) when the interband pairing gab is 
stronger. 
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D. Temperature dependence of and A;, 



Figures 5-7 shows the temperature dependence of and A;,. [27] In Fig. 5, when the 
impurity concentration (uj) increases, \Aa\ and \Afj\ approach each other irrespective oi 9ab- 
Although we do not show the result for < \vab\ < 2 in Fig. 6, the same tendency as in Fig. 
5 is obtained with increasing \vab\ up to \vab\ ^ 1- However, when \vab\ ^ 1, Fig. 6 shows 
that difference between |Aa| and |Aft| becomes more remarkable for larger \vab\ because of 
the same reason as that discussed in Sees. 4.2 and 4.3. This tendency is also obtained in 
Fig. 7 as aj increases. 

Figure 5(b) shows that A^ becomes negative when u becomes strong. This is because the 
Cooper-pair acquires the phase 29ab {— tt) in the tunnelling by interband scatterings: When 
this tunnelling occurs frequently, it is favourable that the phases of the order parameters 
differ by tt in order to avoid the depairing effect caused by the tunnelling of the Cooper-pairs 
having the phase tt. This "7r-phase" disappears when the interband scattering becomes very 
strong: In Fig. 6(b), Af, is positive at Vab = 10. In more detail, the stability of this 7r-phase 
is summarized in Fig. 8: The vr-phase continues to exist when the impurity concentration is 
high {uj > 2) (Fig. 8(a)), while it disappears when the inter- and/or intra-band scatterings 
become strong {\vab\ ^ 3, aj > 1.5), as shown in panels (b) and (c). 

E. Superconducting density of states 

Effects of the impurity concentration and the phase / magnitude of the impurity scattering 
on the mixing of the order parameters have been clarified in Sec. 4.4. In this section, we 
study how the mixing effect is actually observed in the superconducting density of states 
(SDOS), because the tunnelling measurement is useful in observing this phenomenon. 

The superconducting density of states in the j-band is given by 

N,{u;) = N,me[ ^ ] = M^=^]- (4-24) 

^Z,((^)2-A^.(^)2-A^.((^)2 V^-^^- 
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Here Zj{u!) — iu!^{iu!jn — > cu+iS) is the analytic continued renormalized Matsubara frequency 



~ R 



in eq. (4.4); A^- (a;) and A^ (u;) are, respectively, the analytic continued quantities of the 
real and imaginary part of Aj in eq. (4.5). Their expressions are given by (Aj(a;)^ = 



A„ (o^) = Re[AJ + iT^[ r°^L (..)cos(2e„,)-A,(..)sm(2e„,) a„ j .^^25) 



where 



Z„(u;)^6(c^) - (A^A^ + A^AJ cos(2^„,) - (A^AI - A^A J sin(2e„,) 

1 + a, 



2 



^Z^icur - A„(a;) Vz,(a;)2 - A,(a;)2 



(4.26) 



The valuables in the 6-band {Zb{uj), A^ (cu), Af,(ci;) and Db{uj)) are also given by eqs. (4.25) 

and (4.26) where the band indices 'a' and '6' are interchanged. In calculating SDOS, it is 

j ~ R ^7 

convenient to use the last expression in eq. (4.24), where Ej — y Aj (o;)^ + Aj{ujy/Zj{uj) = 



Sf(a;)2 + Sj(a;)2: 

Rf^rA 1 I i "g l^^cbP Sf cos(2e„b)-S^sin(26'„6) 



TrnfA 1 M »a |^,,|2 E:«sm(2e,6)+S^cos(2e,i) 



a; + i- 



\Vab\^ 



Hf and are obtained from eq. (4.27) by interchanging the indices 'a' and '6'. 

Figures 9 and 10 show impurity effects on SDOS in the cases of 9ab — and Oab — T^/'^, 
respectively. When the impurity concentration increases, SDOS's in the two bands are mixed 
with each other to be almost the same at = 20 in Fig. 9(a) and u — b in Fig 10(a). In 
Fig. 10(a), because of the strong depairing effect by the phase of the interband scattering 
9ab^ the mixed SDOS at m = 5 has a small energy gap compared with the case of m = 0.1. 
This strong depairing effect also leads to a gapless behavior at n = 1 in Fig. 10(a). 
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The same mixing effect is observed also in Figs 9(b) and 10(b): Although SDOS's in the 
two bands are different at \vab\ = 0.1, they become almost the same at \vab\ = 1; however, 
since the interband impurity effect again becomes weak when \vab\ ^ 1, SDOS's differ from 
each other at \vab\ — 100. Since the effect of the interband scattering is also weakened by 
the intraband scattering, we obtain two distinct SDOS's at a 10 in Figs. 9(c) and 10(c). 

If MgB2 is really a two-band superconductor with two distinct order parameters, we 
expect the mixing effect in SDOS shown in Figs. 9(a) and 10(a) as the impurity concentration 
is increased. This effect is independent of the phase of the interband scattering as shown 
in Figs. 9 and 10 and, in principle, it always occurs when the impurity concentration is 
high enough. This effect might be useful to judge whether or not superconductivity in 
MgB2 has two distinct order parameters in the two bands. However, infinitely high impurity 
concentration is impossible experimentally. Thus when we try to observe the mixing effect, 
the character of impurity potential is crucial: The intraband scattering must not be strong, 
and the strength of the interband scattering must be moderate {\vab\ ~ 1)- 



V. SUMMARY 



In this paper, we have investigated nonmagnetic impurity effects on two-band s-wave su- 
perconductivity with two different order parameters. We showed that a bound state appears 
around a nonmagnetic impurity depending on the character of the interband scattering po- 
tential. We also considered the case of finite impurity concentration focusing on the mixing 
of the two order parameters by interband scatterings. In contrast to the previous works 
which treat the impurity scattering within the Born approximation, we took into account 
multi-scattering processes using the i-matrix approximation. We showed that the mixing 
of the two order parameters occurs when the impurity concentration becomes high, as dis- 
cussed previously. However, although the interband scattering is essential for the mixing 
effect, we clarified that this effect disappears when the interband scattering is very strong 
{\vab\ 3> 1.). In addition, the intraband scattering also suppresses the mixing effect. When 
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the magnitude of the impurity potential is very large, a clear two-gap structure is observed in 
the superconducting density of states as in the case of the clean system (unless the impurity 
concentration is very high.). 

In order to observe the mixing effect of the order parameters in a two-band superconduc- 
tor, we have to choose impurities with weak intraband scattering and moderate interband 
scattering {\vab\ ^ !•)• If can choose this kind of impurity for MgB2, it is expected 
that, if the observed two-gap structure in the superconducting density of state originates 
from two order parameters in the two bands, the two-gap structure would gradually change 
into a single-gap structure as the impurity concentration increases. The observation of 
this impurity effect may be useful in clarifying the origin of the two-gap structure in the 
superconducting density of states in MgB2. 
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Appendices 



APPENDIX A: EQUATION OF BOUND STATE ENERGY 



From eq. (3.3), the energy of the bound state u) is determined by the equation 

= det[l - VaailOJm) E <^a(, i^m)] 

P 

= det[l - VabPbVhaPa]]^Q^[^ " ^aa ^a(' i^m)], 



(Al) 



where Pj = iu^m)[l — %i Z]p i^m)] ""^ (j = a, Since Vaa is the intraband 

nonmagnetic impurity potential, no pole is obtained from the factor det[l — Hp i<^m)] 
in eq. (Al). Hence we drop this term in what follows. The factor Pj can be calculated using 



EG'; = -7r7V^(0)- 



Qj 



(A2) 



Here Qj — oJ^ + and we have taken the order parameters real as Kj — Ajpi. Then we 
obtain 



7r7V.(0) 



1 + a] 



(A3) 



Substituting eqs.(A3) and (2.8) into eq. (Al), we obtain the equation of the pole as 

\ 

= det =1 + \A\^ + \B\^ - {A + A*), 



\ 



1 

1 



J 



B* I- A* 

where the matrix elements A and h are given by 

A^ 
B = 



(A4) 



^b 



\Vab\ 
\Vab 



2^aA^m , \ , -2 ^b A^m x 
+ Oib)+Vab-Fr{—-aa) 



(A5) 



Qa ^b " ^6 

Calculating each term in eq. (A4) and taking the analytic continuation icu^ cu + iS, we 
obtain eq. (3.5). 
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APPENDIX B: DERIVATION OF THE RENORMALIZED VARIABLES uij AND 



We consider the a-band. The self-energy in eq. (4.2) can be rewritten as 

1 



^Na (0) Vaa TzNa (0) ^« ( ' ^^ru) 

The summation of the Green function in terms of gives 



(Bl) 



1 



ICO" 



7o + 7iPi + 72P2, 



sin(0a)P2 



(B2) 



where cpj is the phase of Aj {j — a,b). 

Using the expression of t4a in eq. (4.2) and eq. (B2), we obtain 



nNa{0)Vaa = 



\Vab\ ^^m 



COs(26la6 + 0b)pi 



1 + al Qb l + al fib 

, \Vab\'^ |Af,| . /r,/) , 7 N , / ttfe I- |2\ 

1 + at \lb i- + ctt 



"■b 

= Ao + Aipi + A2P2 + A3P3 



(B3) 



Thus 



AiPi + A2P2 + A3P3 - Ao 



(B4) 



where 77a = A^ + A| + A| — Aq is reduced to eq. (4.6). Substituting eqs. (B2) and (B4) into 
eq. (Bl), we have 



u„ 



^aai}^m) = JtH^O + Valo) + (Al - 77a7l)Pl + (^2 " ?7a72)P2 + A3P3]. (B5) 

Here = [(Ai - r]a-fif + (A2 - 77a72)^ + ^1 - (Aq + 77a7o)^]/?7a equals eq. (4.7). Thus from 
eqs. (4.1) and (B5), we find 



(B6) 





= i'^m - ^(Ao + 77a7o), 








COS(0„) = |A„| COs(0a) - 




- Vail) 


v 


sin(0a) = |A„|sin(0„) + 
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where cpj is the phase of Aj. Substituting the detailed expressions of 7^ and Aj into eq. (B6), 
we obtain eqs. (4.4) and (4.5). 

We briefly mention that the term proportional to in eq. (B5) gives a renormalization 
of the kinetic energy as 

This renormalization may be crucial when we consider a correlation function; however, since 
it does not affect the integration in terms of the kinetic energy within the treatment in this 
paper, we can safely neglect it. 

APPENDIX C: MIXING OF THE ORDER PARAMETERS WITHIN THE BORN 

APPROXIMATION 



We consider the case of aa — CKb — gab — Qha — 0. [19] In this case, since the phase of 
the interband scattering 6ab is irrelevant, we can take Vah real. In addition, Aj = Aj and Af, 
can be also taken real. When we retain scattering processes within the Born approximation, 
eqs. (4.4) and (4.5) are reduced to 



(CI) 



A, 



r 



A, 

"a' 
' 



where Fa and are defined in eq. (4.22) where we take 
— — 1. Equations (CI) and (C2) can be rewritten as 



(C2) 



= 77a = ^& = and 



i + r< 



(C3) 
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Here we have introduced o)^ = u-^/ujjn-, = ^jl^j and X = A;,/Aa. When we define 
Yj = A^/t<}^, we find a relation Ya — XYb in the dirty hmit (Fa » 1 and ^ 1). Using 
this relation, we find that Ya and Yf, obey the following equation: 

/ 



I 4- La 



\ xa 



1 + 





(y] 


_M 









(C5) 



where Qr = V^^^^T+^^^Af. The solution of eq. (C5) is given by Yj = A^y/A^ (in the dirty 
limit), where 



7Va(0) + 7V,(0) 

is an averaged order parameter. Then the gap equation (4.8) is reduced to 



(C6) 



Substituting the explicit form of Yj into eq. (C7) we obtain 

9aaNb{0) 



(C7) 



/ ^ 9aaNa{0) TT ^ 



TT 



7V„(0) + 7V,(0)/3 



^&6A^a(0) TT 



A2 



E 



1 - 



E 



E 



^™ a/o;^ + A^y 
1 



1^ 7V„(0)+7V,(0)/?t:^a;^ + Aiy A^a(O) + A^.(O) /3 ^(^^ + A^^ ^ 



VA,^ 



Equation (C8) has a solution when the following equation is satisfied: 

Equation (C9) can be interpreted as the gap equation of Aav with the averaged pairing 
interaction (^a„7Va(0) + ^;,6A^6(0))/(7V„(0) + 7V,(0)). 
In the dirty limit, we obtain the following relation: 
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"""^ - """^ (C.n.>0). (CIO) 



Since the superconducting density of states in the j-band is obtained from the analytic 
continuation of LHS in eq. (CIO), we find that the excitation gap in both bands is equal to 
A^y. Namely, in the dirty limit, only one excitation gap Aav appears in the superconducting 
density of states within the Born approximation. 

APPENDIX D: EVALUATION OF EQ. (??) 



When we expand eq. (4.16), we obtain 

- 9cui9bb^J2^M^ ■'^Y.^(^'m)^-9ab9ba^^s{uJm)^ " ^I3^(^m)^- (Dl) 

IjJm a ui'^ b Wm a ui'^ b 

The first and the second terms in RHS in eq. (Dl) include the following complex terms: 

2^S[UJm)[9baj;- ^2 ^~57 +g^«''7r i i ^2 O O/ ^ ^^^^ 

However, from eqs. (4.7) and (4.10) and the definition of Uj, we obtain (1 + al)Da — 
(1 + al)Df, and Na{0)ua — Nb{0)uf,. Substituting these relations into eq. (D2) and noting 
that Qab = 9ba and Clj = Cl'j/\F^\, we find that eq. (D2) is real. Thus the first and the second 
terms in eq. (Dl) is also real. 

The third and the fourth terms in eq. (Dl) can be rewritten as 

/- - _- - V _^(^m)g(^^) r.^ Uh _ Vb Ug _r]a 

[9aa9bb ^a.^.aj^2^2.^ 0:,(o;™)nU<)^^' A(a;^) l],(a;^) ^ i^a«)^^.«)^ 



UaUblVabl"^ 



(1 + al){l + al)Dg{u;m)Db{u;'jnaMnb{^'y 
Equation (D3) is real, so that eq. (4.16) is found real. 



(D3) 
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Figure Captions 

Fig.l: Temperature dependence of the order parameters in the clean system. We put Na{0) — 
Nh{0). This condition is used throughout this paper. In this case, (jab = Qba is satisfied. 
In panels (a) and (b), the intraband interactions gaa and gbb are attractive; panel (a) 
shows the case of gab ~ {gaa, Qbb), while panel (b) is the case of gab < {gaa, Qbb)- Panel 
(c) shows the case when gbb is repulsive. 

Fig. 2: Dependence of the bound state energy on the interband impurity scattering Vab- We 
take Afl = 2, A;, = 1 and = 0. Panel (a) shows the effect of the phase of Vab in the 
case of ab — 0. Since A^ — cos{29ab)^a^b > for 9ab > 7r/6, the bound state is absent 
for 9ab > 7r/6. Panel (b) shows the effect of the intraband potential scattering ab in 
the 6-band in the case of 6'a(, = 7r/2. 

Fig.3: Dependence of the transition temperature on the interband potential scattering Vab- 
The transition temperature is normalized by the value at Vab — 0. We take gaa — 0.3, 
^66 = 0.25, Ua — Ub — 1 and = 0. (a) Effect of the interband pairing interaction gab 
in the case of 6ab — and Ub = 0. (b) Effect of the phase of the interband potential 
scattering 6ab in the case of gab — 0.002 and ab = 0. (b) Effect of the intraband 
potential scattering in the 6-band {ab). We put gab — 0.002 and 9ab — 0. 

Fig. 4: Dependence of the transition temperature Tc on the impurity concentration Uj — 
nimp/{7^Nj{0)). In this figure, we take — Ub and is normalized by the value 
in clean system (T°). The other parameters are taken as gaa = 0.3, gbb = 0.25, 
Oia — Oib — Q and \vab\ — 1- In panel (a), the phase of the interband scattering is taken 
6ab — 0, while in panel (b) 6ab = 7r/2. 

Fig. 5: Temperature dependence of A^ and A^ in the cases of (a) 9ab — and (b) 9ab — 7r/2. 
We take gaa = 0.3, gbb = 0.25, gab = 0.002, aa = ab = and \vab\ = 1- The x- and 
y-axis are normalized by Tc in the clean system. 
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Fig.6: Effect of the interband scattering on temperature dependence of and A;, in the 

cases of (a) 9ab = and (b) 9ab = 7r/2. We take Qaa = 0.3, gbb = 0.25, gab = 0.002, 
Q^a = ctfe = and Ua = Ub = 1. 

Fig. 7: Effect of intraband scattering («„ and a^) on T^. The figure shows the case of Oab = 
{vab — 1)- The other parameters are Qaa — 0.3, gbb — 0.25, gab — 0.002, Ua — Ub — 10 
and \voi\ = 1- 

Fig.8: Stabihty of "7r-phase" caused by the interband scattering with Oab — '^1'^- We take 
T = 0.5Tc, gaa = 0.3, gbb = 0.25 and gab = 0.1. 

Fig.9: Superconducting density of states at T = O.lTc in the case of 9ab = 0. We take 
9aa = 0.3, gbb = 0.25 and gab = 0.002. Nj^uj) is normalized by the density of states at 
the Fermi level in the normal state (Na(0) = Nb(0) = N(0)). In this figure, 'a' and 'b' 
mean SDOS in the a- and 6-band, respectively, (a) Effect of impurity concentration 
{u = Ua = Ub) at = ctft = and \vab\ = 1- Since SDOS's are almost the same in 
the two bands a,t u — 5, we show one of them only in the figure, (b) Effect of the 
magnitude of the interband scattering \vab\ at = a?, = and Ua — Ub = 10. (b) 
Effect of the intraband scattering a {= ota = otb) at \vab\ = 1 and Ua — Ub= 10. 

Fig. 10: Superconducting density of states at T = O.lTc in the case of Oab = tt/2. We take 
gaa = 0.3, gbb = 0.25 and gab = 0.1. (a) Effect of u Ua ^ Ub) at Q;a = Q;b = 
and \vab\ — 1- (b) Effect of \vab\ at cca = = and Ua — Ub — 5. (c) Effect of 
the intraband scattering a at \vab\ = 1 and Ua = Ub = 5. Since Na{uj) and Ni,{u) are 
almost the same at u = 5 in panel (a), \vab\ — 100 in panel (b) and a = in panel (c), 
only Na{u;) is shown in these cases. 
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Fig. 5 Y.Ohashi 
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Fig. 6 Y.Ohashi 




Fig.7 Y.Ohashi 
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Fig.8 Y.Ohashi 
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